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ABSTRACT 

The density irregularities and holes visible in many Chandra X-ray images of cluster and galactic 
cooling flows can be produced by symmetrically heated gas near the central galactic black hole. As the 
heated gas rises away from the galactic center, a relatively small number of large plumes and bubbles are 
formed in qualitative agreement with the observed features. The expanding centrally heated gas drives 
a shock into the surrounding gas, displacing it radially. Both computational and analytic results show 
that the ambient gas near the bubble is cooled by expansion, accounting for the cool rims commonly 
observed around X-ray holes in cooling flows. 

Subject headings: galaxies: elliptical and lenticular, CD - galaxies: active - cooling flows - X-rays: 
galaxies - galaxies: clusters: general - X-rays: galaxies: clusters 



1. INTRODUCTION 

High resolution X-ray images of "cooling flows" in ellip- 
tical galaxies taken with the Chandra observatory indicate 
that the central hot gas in these systems is not smoothly 
distributed, but is cavitated on scales comparable to the 
radio emission (M87: Bohringer et al. 1995; Fabian et al 
2001; Wilson et al 2001; M84/NGC 4374: Finoguenov & 
Jones 2001: NGC 4636: Jones et al. 2001; Loewenstein 
et al. 2001; David et al. 2002; NGC 507: Forman et al. 
2001; NGC 5044: Buote et al. 2002). In view of their 
short dynamical times, 10 7 — 10 8 yrs, these disturbances, 
and by inference also the radio sources, must be highly 
transient. X-ray holes are also seen in galaxy clusters on 
scales of < 50 kpc, often approximately coincident with 
lobes of extended radio emission (e.g. Hydra A/3C295: 
McNamara et al. 2000; David et al. 2001; Allen et al 
2001; Perseus/NGC 1275/3C84: Churazov et al. 2000; 
Bohringer et al. 1993; Fabian et al 2000; A2052: Blanton 
et al 2001; A4059: Heinz et al 2002). The Perseus cluster 
in particular contains a multitude of X-ray holes located 
at random azimuthal orientations. The X-ray holes are 
regions in which the gas is either heated or displaced by 
a non-thermal plasma of comparable energy density. In 
either case the low-density holes must be buoyant (e.g. 
Churazov et al. 2001). 

One of the most remarkable features of these holes is 
that the X-ray gas around the rims is often colder than the 
average nearby ambient gas (e.g. Fabian et al. 2000, 2001; 
Fabian 2002). Low temperature rims have been held as 
evidence that the surrounding gas has not been (recently) 
strongly shocked, although shocks must accompany even 
a subsonic expansion of the hot gas. These low entropy 
rims cannot be understood as local gas that was shocked 
and subsequently lost entropy by radiation (Nulsen et al 
2002; Soker, Blanton & Sarazin 2002), but may instead be 
low-entropy gas that has somehow been raised from the 
center of the flow as these authors suggest. 



The random location of older, more distant bubbles in 
Perseus and the irregular X-ray images within the central 
few kpcs of elliptical galaxies, are apparently inconsistent 
with non-thermal jets having fixed orientations defined by 
the spin axis of massive black holes. Conversely, not ev- 
ery strong radio lobe corresponds to an X-ray hole, as in 
Perseus/NGC1275/3C84. 

We show here that both the random orientation of X-ray 
holes and their cold rims can be understood if low entropy 
hot gas near the center of a cooling flow is heated by an 
active nucleus in the central elliptical galaxy. 

2. PREVIOUS STUDIES OF X-RAY HOLES 

Several recent theoretical studies of X-ray bubbles have 
examined the consequences of introducing heated gas in 
some localized region away from the center of the cool- 
ing flow. The 2D hot buoyant bubbles by Churazov et 
al. (2001), for example, are accompanied by a column of 
colder (low entropy) gas that moves radially upward near 
the center of the bubble. These bubbles float upward in 
the cooling flow atmosphere and come to rest at some large 
radius where the ambient entropy is equal to that within 
the bubble. According to Churazov et al, rising bubbles 
capture cold gas and lift it to large distances from the 
center, although if this gas is not heated it would be ex- 
pected to eventually fall back. Nevertheless, the heating 
required to form the bubbles reduces the global entropy 
gradient in the cooling flow. Quilis et al. (2001) studied 
the evolution of a nearly axisymmetric 3D bubble pro- 
duced by heating gas at some finite radius. They noticed 
that the buoyant bubble is surrounded by a shell of slightly 
colder gas which they attributed to cooling expansion as 
the gas is pushed by the bubble toward regions of decreas- 
ing ambient pressure. In these studies gas is assumed to be 
heated at some off-center site, but the implication is that 
the ultimate source of energy is located in a central active 
black hole. Heating by jets that proceed directly from the 
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AGN core can also create buoyant regions and upwelling 
of low entropy gas, but only along the jet axis (Reynolds, 
Heinz, & Begelman 2001; 2002). Reynolds et al. (2002) 
also note that "curiously, the [observed] cool gas seems to 
form a shell around the radio lobes and occupies precisely 
the location where we would expect shock/compressionally 
heated gas to reside" . Regardless of the heating geometry 
assumed, by jets or ad hoc off-center bubbles, the global, 
long-term rate that hot gas cools by radiative losses in 
cooling flows is unlikely to change without altering the 
gas density and temperature profiles in disagreement with 
observations. (Brighenti & Mathews 2002). 

3. PRODUCING BUBBLES WITH CENTRAL HEATING 

We begin by constructing a simple cooling flow for a 
massive elliptical galaxy. The gravitational potential is a 
superposition of a de Vaucouleurs stellar mass profile of 
total mass M* t = 7.26 x 10 11 M and a Navarro- Frenk- 
White dark halo with virial mass M h ^ vir = 4 x 10 13 M 
and concentration c = 10, typical for the galaxy group 
environments in which massive E galaxies are thought to 
form. The initial cooling flow that we consider is gener- 
ated entirely by mass loss from an old stellar population, 
a*p* = 4.7 x 10 -20 (t/t n ) -1 - 3 gm s _1 where p*(r) is the 
stellar density and t n = 13 Gyrs is the assumed age of the 
stars. Some heating is provided by Type la supernovae, ex- 
ploding at a rate of 0.06(t/t n ) -1 per 100 years per 10 10 L#, 
each releasing 10 51 ergs. Supernova heating at this level is 
consistent with the iron abundance observed in the hot gas 
but has a relatively minor influence on the overall energet- 
ics of the hot gas. The gas-dynamical equations that we 
solve are identical to those described in Brighenti & Math- 
ews (2002). We use a 2D code modeled after ZEUS (Stone 
& Norman 1992) with 490 grid zones in both cylindrical 
coordinates. For R, z < 50 kpc the zones are uniform with 
size 125 x 125 pc. Outside this region the size increases 
geometrically toward the outer boundary at 260 x 260 kpc. 

The quiescent cooling flow that results after slowly 
evolving for 13 Gyrs has a characteristic temperature of 
~ 1 keV and a gas density profile p oc r -1,2 within about 
10 kpc and slightly steeper beyond. We assume that gas 
within rh = 1 kpc of the center is heated at a constant 
rate L h = 2.5 x 10 42 erg s" 1 for t h = 5 x 10 7 yrs. The 
heated gas reaches temperatures far in excess of 1 keV 
and has a sound crossing time much less than th, allowing 
an ultra-hot low-density region to develop near the galac- 
tic core. The hot gas pushes outward, driving a shock 
into the cooling flow gas beyond. The resulting evolu- 
tion of the gas density and temperature are shown in Fig- 
ure 1 at three times t x = 40 x 10 6 , t 2 = 80 x 10 6 , and 
t 3 = 100 x 10 6 yrs. The Ray leigh- Taylor (RT) develop- 
ment of a few large plumes filled with heated gas produces 
irregular X-ray patterns that resemble Chandra images of 
NGC 5044 and other E galaxies. Small wavelength RT un- 
stable features are not prominent even though they would 
be resolved by the numerical grid. The hot cavity expands 
gently and is already subsonic (< 500 km s _1 ) at t = 1 
Myr. By time t\ the weakened precursor shock has moved 
out to ~ 20 kpc, moving almost sonically. 

Low temperature, bright rims surrounding the hot 
plumes/bubbles are clearly visible in Figure 1. The peak 
density in the rims at t\ is several times that of the ini- 



tial unshocked gas at this radius. The minimum entropy 
log^/n 2 / 3 ) in the simulation at times t = 0, ti, t 2 and 
ts are 7.15, 7.9, 8.1 and 8.0 respectively. Low entropy gas 
in the rims at t > t\ was originally at r ~ as verified 
with passively advected tracer particles. This gas cools by 
adiabatic expansion while rising in the hot gas atmosphere 
(see also Nulsen et al. 2002; Soker et al. 2002). The larger 
bubble at time t 2 is penetrated by a disorganized S-shaped 
column of slowly rising cold gas. The upper panel of Fig- 
ure 2 shows the X-ray surface brightness of the central 
cooling flow at time t 2 when the plumes have separated 
into clearly visible buoyant X-ray bubbles. 

It is important to stress that the physical size of the 
plumes that form is not determined by our computational 
resolution. Calculations repeated at higher spatial resolu- 
tion (central grid size=50 pc) also produce a few plumes 
and (ultimately) bubbles of about the same size as those 
shown in Figures 1 and 2, which are similar to the X-ray 
patterns observed. An examination of the effective gravity 
in the gas rest frame indicates that the outer rims of the 
large bubbles are RT unstable. At time ts the gravita- 
tional freefall time across the bubble radius is tff w O.U3; 
coherent hot bubbles persist for times ^ tff, even at our 
highest computational resolution. 

We performed similar calculations on more extended 
cooling flows in galaxy groups and clusters with very sim- 
ilar results. For example, the lower panel of Figure 2 illus- 
trates the X-ray surface brightness at time t = 120 x 10 6 
years produced in a rich cluster cooling flow (M vir = 10 15 
Mq and T vir ^6.9 keV) with heating parameters: Lh = 
10 45 ergs s _1 , rh = 2 kpc and th = 5 x 10 T yrs. Very little 
gas cools by radiation losses during the brief duration of 
our calculations shown in Figures 1 and 2 so the flow is 
essentially adiabatic. This is expected since a defining at- 
tribute of cooling flows is that the dynamical time is much 
less than the radiative cooling time. 

4. SIMILARITY SOLUTION 

Cold rims are a generic feature of flows forced to ex- 
pand into a medium of increasing entropy. This can be 
illustrated with self-similar flows described by the vari- 
able n = r~ x t (Courant & Friedrichs 1948; Rogers, 1957; 
Parker, 1961; Chevalier & Imamura 1983). The adiabatic 
equations for spherical flow can be written in terms of 
r] provided u = (r/t)U(rj), p = A r~ a Vt(j]) and P = 
A r 2 ~ a t~ 2 U(r]). The hot region of centrally heated gas 
can be represented by an expanding spherical piston at 
r] = r] 2 that drives an advancing shock at n = rji < n 2 . The 
solutions U(rj), Q(rj) and 11(77) can be found by inward in- 
tegration from post-shock values J7i, Oi and II 1 which are 
functions of A, 7 = 5/3 and the Mach number of the shock 
M (e.g. Rogers 1957). We consider the spherical flow 
that results for pistons uniformly expanding (A = 1) into 
an isothermal, fully ionized medium (T Q = 10 T K) with 
decreasing electron density n e = 0.0526r _Q! cm -3 where 
a = 1.2 approximates the gas density observed in mas- 
sive elliptical galaxies (e.g. Brighenti & Mathews 1997). 
Figure 3 illustrates the flow generated by three pistons 
driving shocks at r\ = 3 kpc with M. = 1.5, 2 and 3. In 
each case the temperature drops sharply near the piston, to 
values even less than T , creating cold rims where the den- 
sity rises accordingly. Post-shock gas that has just cooled 
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back to T Q in Figure 3, located at \og(r/r\) = —0.179, 
-0.108, and -0.078 for M = 1.5, 2 and 3 respectively, 
was originally located at log(r/ri) = —0.347, -0.476 and 
-0.957, much closer to the origin. At the piston where 
U < 1 approaches unity, the gas temperature decreases as 
X = jU/Q oc (1 - Uy where p = 2a/3(5 - a) = 0.21 when 
A = 1. The cold gas adjacent to the piston is a relic of low 
entropy gas initially near the center that was shocked then 
cooled by expansion; d\og p/d\og t is negative throughout 
the flowing gas in the plotted solutions. Similarity flows 
with lower shock Mach numbers have flatter density pro- 
files; for M = 1.5 there is a shallow density minimum at 
log(r/ri) = —0.047, indicating weak post-shock RT insta- 
bility. However, the constant piston velocity assumed in 
the similarity flows differs from the rapidly decelerating 
heated bubble in the computed galactic scale flow. 

5. FINAL REMARKS AND CONCLUSIONS 

If central heating is responsible for the frequently ob- 
served X-ray holes and surface brightness fluctuations in 
typical Chandra images, as we claim, then such AGN-black 
hole heating may be a nearly universal component of cool- 
ing flows. The physical origin and nature of this heating 
- e.g. shock waves, cosmic rays, etc. - has not yet been 
explored but these details are not essential to the forma- 
tion of plumes and bubbles. While the production of X- 
ray irregularities is quite generally insensitive to the as- 
sumed heating parameters, we have noticed several trends 
when Lh, rh and th are varied. Non-linear plumes mature 
faster when is reduced. If r^ is too large, the cold rims 
are less pronounced. However, the possibility of learn- 
ing about such heating details by comparing X-ray images 
with gasdynamical calculations is limited if the mean time 
between heating episodes is sufficiently short that the ini- 
tial pre-heated gas is no longer perfectly quiescent as we 
have supposed here. 

All well-observed cooling flows emit diffuse optical 
Balmer and forbidden lines visible in the brightest central 
regions. Enhanced Ha + [Nil] emission has been observed 
near the boundaries of X-ray bubbles (e.g. Blanton et al 
2001; Trinchieri & Goudfrooij 2002). It is tempting to 
conclude that such gas has cooled from the hot phase in 
the cooler, denser gas around the rims of the plumes or 
bubbles. But direct cooling from the hot phase is not sup- 
ported by our detailed hydrodynamic calculations which 
are essentially unchanged if radiative cooling losses are ne- 
glected. The radiative cooling time is greater than the to- 



tal time of our simulations. The cooler gas at T ~ 10 4 K 
responsible for the observed optical line emission may sim- 
ply be displaced by the plumes and bubbles along with the 
local hot gas. However, on longer timescales than we con- 
sider here, heating in the central regions stimulates cooling 
by thermal instabilities in more distant regions of the flow, 
so the time-averaged global cooling rate is essentially un- 
changed (Brighenti & Mathews 2002). 

If intermittent AGN heating is reasonably frequent, the 
entire inner cooling flow will become convective and tur- 
bulent (Brighenti & Mathews 2002). Turbulent motions 
observed in the diffuse optical lines (e.g. Caon et al. 
2000) are too energetic on large scales to be understood 
in terms of supernova explosions, but could result from 
central heating similar to that described here. 

X-ray holes formed by central heating are randomly dis- 
posed around the center of the flow, as in the Perseus 
cluster. This is a particularly desirable feature since the 
X-ray bubbles need not be aligned with the accretion axis 
of the central black hole. Nevertheless, radio lobes often 
appear to be associated with X-ray holes. This connection 
may be causal, but in some cases the relativistic electrons 
may readily flow into pre-existing holes; this may be oc- 
curring in the central southwestern hole in Perseus (Fig. 
1 of Fabian 2001). 

Finally, the azimuthally averaged hot gas density pro- 
file in these centrally heated models is flattened by the 
holes, providing an excellent fit to observed X-ray surface 
brightness profiles which are too centrally peaked other- 
wise. To maintain this fit, AGN heating would need to 
reoccur every ~ 10 8 years. However, in the presence of 
central heating the gas temperature profile is less satisfac- 
tory, lacking the deep central minimum usually observed. 

We conclude that central heating in cooling flows natu- 
rally results in a relatively small number of randomly lo- 
cated X-ray bubbles similar to the irregularities observed. 
In addition, these expanding plumes and bubbles are sur- 
rounded by considerably cooler gas, also consistent with 
recent Chandra observations. Gas in the cold rims was 
shock-heated but subsequently cooled by expansion. 
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Fig. 1. — Evolution of centrally heated gas in a galactic scale cooling flow. The z (vertical) and R (horizontal) coordinates in each panel 
are labeled in kpc. The upper and lower series of panels show respectively the gas density and temperature at three times: t\ = 40 x 10 6 , 
£2 = 80 x 10 6 , and £3 = 100 x 10 6 yrs from left to right. Already at time £1 the initial precursor shock has moved out to ~ 20 kpc and 
additional shocks are seen at time £3. The hot plumes and bubbles are surrounded by relatively cold gas, particularly at early times. This 
figure can be viewed in color at the electronic ApJ site. 

Fig. 2. — X-ray surface brightness distributions (0.5 - 2 keV) showing bubbles in two flows. The z (vertical) and R (horizontal) coordinates 
in each panel are labeled in kpc. Upper panel: For the galactic scale cooling flow in Figure 1 shown at time £2 = 80 x 10 6 yrs. Lower panel: 
For a rich cluster cooling flow at time £ = 120 x 10 6 yrs with heating parameters: = 10 45 ergs s _1 , r^ = 2 kpc and = 5 x 10 7 yrs. This 
figure can be viewed in color at the electronic ApJ site. 

Fig. 3. — Self-similar flow of gas produced by a spherical piston expanding uniformly into an isothermal, fully ionized static gas with initial 
density n e = 0.0526r — a cm -3 . In descending order the panels show the radial dependence of the electron density, temperature, pressure and 
velocity. The flow is shown at the (arbitrary) moment when the shock has arrived at r\ = 3 kpc with Mach numbers 1.5, 2 and 3. The radial 
variation of the initially undisturbed gas is shown with dashed lines. As suggested by the plots, the temperature in the flowing gas near the 
piston (at the left of each plot) decreases to zero as the density becomes infinite. 
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